Response to tilted magnetic fields in Bi 2 Sr2CaCu20s with columnar defects 
Evidences for transverse Meissner effect 
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The transverse Meissner effect (TME) in the highly layered superconductor Bi2Sr2CaCu20s+ y 
with columnar defects is investigated by transport measurements. We present detailed evidence 
for the persistence of the Bose glass phase for H± < H± c : (i) the variable-range vortex hopping 
process for low currents crosses over to the half-loops regime for high currents; (ii) in both regimes 
near H± c the energy barriers vanish linearly with H± ; (iii) the transition temperature is governed 
by T B a(H h 0) - T B a(H\\,H_i) ~ \H ± \ 1/vj - with v x = 1.0 ± 0.1. Furthermore, above the transition 
as H± — > H± c , moving kink chains consistent with a commensurate-incommensurate transition 
scenario are observed. These results thereby clearly show the existence of the TME for H± < H± c . 
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The interplay between elasticity, interactions, thermal 
fluctuations and dimensionality of magnetic flux lines 
subject to pinning, yields a large variety of vortex phases 
in high-T c superconductors [Q. The nature of these 
phases and how they depend on the pinning is far from 
being completely elucidated. The theory of pinning by 
correlated disorder, such as twin planes or amorphous 
columnar tracks created by heavy ions, has been consid- 
ered by Nelson and Vinokur (NV) [^J and Hwa et al ||. 
In the case of parallel tracks, NV have shown that if the 
applied magnetic field H is aligned with the columnar 
defects, the low temperature physics of vortices is simi- 
lar to that of the Bose glass (BG) j|, with the flux lines 
strongly localized in the tracks leading to zero creep re- 
sistivity. When H is tilted at an angle 8 away from the 
column direction, the BG phase with perfect alignment 
of the internal flux density B parallel to the columns is 
predicted to be stable up to a critical transverse field 
H± c producing the so-called transverse Meissner effect 
(TME) H §. For 8 > 8 C = arctan(^ # ic /B), the lin- 
ear resistivity is finite resulting from the appearance of 
kink chains along the transverse direction as discussed in 
Ref. H Finally, above a still larger angle, the kinked vor- 
tex structures disappear and B becomes collinear with 
H. Similar scenarios apply to vortex pinning by the twin 
boundaries p ] and by the layered structure of the coum- 
pound itself [A, 

Recently, the TME in untwinned single crystals of 
Y Ba 2 Cu 3l O'j (YBCO) with columnar defects, has been 
observed using Hall sensors Previous magnetic mea- 
surements in anisotropic 3D superconductors, such as 
YBCO, gave support to the presence of vortex lock- 
in phenomena, due to pinning by the twin boundaries 
H or by the interlayers between the Cu-0 planes |jio| . 
In the case of highly layered superconductors, such as 
Bi2Sr2CaCu20s+ y (BSCCO), the lock-in transition was 
observed for H tilted away from the layers jll), but the 
existence of the TME due to the pinning by columnar 
defects remained an open question. 

In this letter, we present measurements of the electrical 



properties near the glass-liquid transition in BSCCO sin- 
gle crystals with parallel columnar defects, as a function 
of T and 8 for filling factors / = [iqHu/Bq, < 1. Below 
a critical angle 9 C (T), we observe a vanishing resistivity 
p(J) for low currents. A detailed quantitative analysis of 
p{ J) indicates that the creep proceeds via variable-range 
vortex hopping (VRH) at low currents due to some dis- 
order iQ, crossing over to the half-loop (HL) regime at 
high currents. For 8 > 8 C (T), the very signature of a 
kinked vortex structure, consistent with a commensurate- 
incommensurate (CIC) transition scenario in (1+1) di- 
mensions H , is deduced from the critical behavior of the 
linear resistivity. All these results clearly demonstrate 
that when H is tilted at 8 < 8 C away from the defects, 
the flux lines remain localized on columnar defects, and 
hence, the BG exhibits a TME. 

The BSSCO single crystal was grown by a self-flux 
technique, as described elsewhere fl3fl . The crystal of 
1 x 1 x 0.030 mm 3 size with the c-axis along the shortest 
dimension has a T c of 89 K, a transition width of ~ 1 K, 
and was irradiated along its c-axis with 5.8 GeV Pb ions 
to a dose corresponding to _B$ = 1.5 T at the GANIL 
(France). Isothermal I — V curves were recorded using 
a dc four-probe method with a sensitivity of ~ 10~ 10 V 
and a temperature stability better than 5 mK . H was 
aligned with the tracks using the well-known dip feature 
occurring in dissipation process for 8 = 0°, and was tilted 
with an angular resolution better than 0.1° away from the 
column direction for / fixed. 

Fig.?? shows a typical log-log plot of V/I — I curves 
obtained varying 8 for / = 1/3 and T fixed below the 
BG transition temperature determined at 8 = 0° [0. 
We observe a well-defined angular crossover at an angle 
8 C . For 8 < 9 C , data suggest that the BG phase persists. 
On the contrary, the existence of an Ohmic regime for 
8 > 8 C indicates a liquid vortex-like state. 

We first focus on the angular range 8 < 8 C , and we con- 
sider the scenario where the BG phase is stable. Thus, 
one expects that excitations of some localized vortex lines 
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lead to a nonlinear resistivity given by : 

p{J) = po exp( - E K {J c /jy/k B T) (1) 

where p is a characteristic flux-flow resistivity and 
Ek{Jc/ J)^ represents the barriers against vortex motion. 
This expression is predicted to hold for various regimes of 
different behavior as the current probes different length 
scales in the BG. Ek acts as a scaling parameter for 
the pinning energy and J c is the characteristic current 
scale of the creep process. When the current is large 
enough that the growth of vortex-loops excitations of a 
line from its pinning track does not reach out the neigh- 
boring tracks, the HL excitations are relevant and yield 
an exponent p — 1 and J c = J\ = U/^od), where U is 
the mean pinning potential and d = y^^>o/ i?$ is the mean 
spacing between pins. With decreasing current, size of 
half-loops increases with the result that some disorder in 
the pinning potential becomes relevant. This situation 
yields the VRH process characterized by p = 1/3 and by 
another important current scale J c = Jq = l/( < f > o.g(/t)d 3 ), 
where g(p) denotes the density of pinning energies at the 
chemical potential of the vortex system. One expects 
that both Ji and Jo are insensitive to H±. We therefore 
consider that the only important effect of H± is to lower 
Ek, according to the formula : 



Ek = Ek — e 2 <$> Q dH 1 



where Ek is the mean kink energy, and the energy 
gain due to the tilt is obtained from the isotropic result 
$>odH± by applying the scaling rule [Eq. (3.12)] of the 
Ref. 0. 

We now present the following two-step method of 
fitting Eqs. (1) and (2) to the data, which allows to get 
a good understanding of the physics in our experiment. 
First, we plot in Fig.|] the natural logarithm of V/I 
versus tan# for / fixed. Note that tan (9 is here directly 
proportional to H± since / is kept constant. A linear 
variation ln(V/I) — M + Ntan9 (expression designated 
below as El) with M and N being current-dependent 
parameters is found for 9 varying up to 9 C (solid lines), 
whereas above 6 C data deviate from this behavior (dotted 
lines). This finding is another argument supporting a 
true angular transition at 9 C , as suggested in Fig.??. It 
should be noted that we draw such a conclusion from 
the observation of two independent regimes of different 
behavior. One leads to a vanishing linear resistivity 
as 9 — > and the other is evidenced by probing the 
regions of nonlinear resistivity above and below 9 C . 
Another result is that N ~ I~^ with p undergoing a 
jump from 1/3 to 1 at li 60 mi. Such a current 
crossover is clearly visible in insert a of Fig.|| (filled 
symbols) where we plot together N normalized by 

1/3 

I and N normalized by I\ versus / depending upon 
whether / is < I2 or > I2, respectively. Here, Iq and 
I\ are two currents evaluated on the basis of the BG 



model, as will be seen below. It therefrom follows that 
N 1/3 {I) = K(I /I) 1/3 and Ni(I) = KI X /I where the 
labels 1/3 and 1 differentiate between two fits below 
and above I2, respectively. In both these equations, 
K ps 0.6 is a dimensionless constant. It should be noted 
that such a crossover at I2 is clearly observed from 
V/I — I curves as well (see Fig.??), where a sudden 
increase of V / 1 suggesting an increasing vortex motion, 
occurs at 12- Second, we plot in Fig .0 and its insert the 
natural logarithm of V/I versus I - " 3 for I < I2 and 
versus I^ 1 for I > I2, respectively. We verify that the 
expression ln(V/I) = P — Q^I~^ (expression designated 
below as E2) fits very well our experimental data (solid 
lines) with an exponent p — 1/3 for I < I2 and with 
p = 1 for I > I2. In E2, is labeled using the above 
convention. The result is that P is a constant (within 
the experimental errors), while depends on 9. In 
Fig.|| and its lower insert, we show Q1/3 and Q\ versus 
tan (9, respectively. A linear variation in tan# is observed 
consistently with expression El. In upper insert of 
Fig.||, we plot together Qi/3/I^ 3 and Qi/h versus 
tan 9. The result is that the data are superimposable 
onto a single straight line with slope w —0.6 consistent 
with the value of — K, as may be verified by identifying 
the expression El with the expression El. 

In conclusion, V/I = 



(2) R exp 



is clearly observed 



{K'-K tan 6) (I /I) 1/3 
in our experiment for I < I2 (solid lines in Fig.??) 
while V/I = i? exp [- (K' - A' tan (9) (Ji/J)] is valid for 
data above I 2 (dashed lines in Fig.??). In both fitting 
expressions K and K' depend on / and T. In the case 
of data presented in Fig.??, K w 0.6 and K' w 0.4; 
furthermore Iq = 15 A and I\ — 0.4 A are evaluated 
from BG model (see below), and i?o ~ 1-1 pQ is five 
orders of magnitude lower than the normal resistance. 
Therefore, for 9 < 9 C we have rather strong evidence 
of two separate (albeit related) vortex creep processes 
peculiar to a BG : the HL expansion with p ~ 1 which is 
cut off by the crossover at I2 into the VRH process with 
p = 1/3. Moreover, we note that for each filling factor 
investigated in our experiment (/ = 2/15, 1/3, 2/3), we 
observe a continuation of the VRH process evidenced 
at 9 = 0° 1 12 . We therefore argue that the BG phase 
remains stable up to a critical tilting angle 9 C , as 
predicted by NV §. 

To test the accuracy of the above view, we estimate 
from the theory |], g], first the characteristic current 
scales, secondly the current crossover, and then the 
energy barriers in Eq. (1). J only depends on g(p), the 
density of pinning energies evaluated at the chemical 
potential of the vortex system. Although a form of 
g(x) is not yet available, an estimate of g(p) can be 
done in terms of the bandwidth of pinning energies 7, 
due to the disorder, and hence g(p) w 1/ (d 2 ^) so that 
J « j/((p d). Following Blatter et al Q, 7 = id + 7* 
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where t d w U / y/E K /k B T exp {-y/2E K /k B T) estimates 
the dispersion in the pinning energies resulting from dis- 
order in the position of tracks and 7, arises from some 
on-site randomness. We shall assume random defect 
radii as the main source of on-site disorder. A reasonable 
estimate of 7^ is then given by the width of the distri- 
bution of pinning energies P(Uk) = P{ck)<1ck / 'dUx 
where Uk is the binding energy of a defect with 
radius ck and P{ck) is the probability distribution 
of the defect radii. A realistic pdf is a normal law 
centered at cq — 45 A with a standard deviation of 
6 A, as shown in figure 1 of Ref. [l2| Thus, using 
the formula U K = {4%/ '8tt fi a \ 2 ab )ln[l + {c K /^ ah ) 2 } 
jl], H where X a b and £ a & are respectively the planar 
penetration depth and the planar coherence length, 
we determine Uq =< Uk > representing the energy 
scale for the vortex pinning. Another important energy 
scale corresponding to vortex positional fluctuations is 
T* = (k B co/4£ Gi)y/ln(ao/£ ab ) (T e — T) w here G l is the 
2D Ginzburg number and ao — yj 4>a/B is the vortex- 
lattice constant. Uq and T* completely characterize 
the system of vortices and columnar defects. Thus, it 
can be inferred E K = (d/v^ab) T*y/f(T/T*) where 
f(x) = x 2 /2 exp{—2x 2 ) accounts for entropic effects, 
and U = Uq f(T/T*). Using appropriate parameters for 
BSCCO (A w 18501, £ « 20 A and G t » 0.2) we esti- 
mate Jo ~ 10 9 A/m 2 and 3 X w 2.6 10 7 A/m 2 for / = 1/3 
and T = 68 K. The balance between the barrier energy 
for the VRH process Uvrh = Ek(Jq/ J) 1 ^ 3 and the bar- 
rier for the HL regime Uhl = Ek Ji/J determines the 
crossover current J 2 = (Jx/Jo) l ^Ji = U 3 ! 2 /((j> dy/j). 
Thus, we have J 2 <C J\ <C Jo so that the crossover 
between the VRH process and the HL regime takes 
place at J2. Assuming uniform currents into the sample, 
we find I2 ~ 60 mA in excellent agreement with the 
experiment (Fig.??). A more complete quantitative 
evaluation of the consistency of the data with the model 
is obtained for different values of / and T, as well as 
for different values of d (insert in Fig.||). In particular, 
we find that the effect of d on I2 is in good agreement 
with the prediction I2 cx \/d. Finally, evaluating Uvrh 
and Uhl with Ek from Eq. (2), we obtain using again 
the above usual parameters of BSCCO with e w 1/200 
that E K /k B T w 0.3 and e 2 (/)odH ± /k B T w 0.7 tan for 
/ = 1/3 and T = 68 K, which correspond to the values 
of our fitting parameters K' and K (see upper insert in 
Fig.||). Thus, the experiment is also in good agreement 
with the barriers Uvrh and Uhl theoretically predicted 
below and above J2, respectively. 

Finally, we pass on to the linear regime observed 
above 8 C (see Fig. 1). For 8 > 8 C , Hwa et al [0] predict 
on the basis of the CIC transition, the appearance of 
free moving chains of kinks oriented in the H± direction 
leading to a critical behavior of the linear resistiv- 
ity p ~ (tan# — tan0 c ) 1J characterized by an exponent 
v = 1/2 [v = 3/2] in (1 + 1) [(2+1)] dimensions. This type 



of behavior with v = 1/2 is evident in Fig.|| which shows 
a plot of R/R c versus tan (9 — tan(9 c for different values 
of / (or equivalently H\\) and T < T B g(H\\,0). Here, R c 
is a scaling parameter comparable to i?0: and for / fixed, 
9 C (T) is determined fitting R = R c (t&n8 - tan^) 1 / 2 
to the data. In Fig.^J, the arrow indicates 9 C ob- 
tained in this way. The insert of the figure 5 shows 
t(H h H ± ) = (T BG (H h 0) - T BG (H h H x )/T Ba (Hp0) 
as a function of H± where T B g(H\\,H±_) is obtained 
through the inversion of 9 C and T. It turns out that 
t(H\\,H±) is independent of H\\. The solid line shown in 
insert is a fit to the data following t(Hu,H±) ~ IH^] 1 ^ 1 - 
with v± = 1.0 + 0.1, as recently suggested from nu- 
merical simulations H] and observed in (K, Ba)BiOz 
|l5|. Moreover, such a value of v±_ is excellently consistent 
with the result previously found using the scaling theory 
of the BG transition at Hj_ = [jjj . 

In summary, we demonstrate from transport mea- 
surements the stability of the BG phase in BSCCO 
single crystal with columnar defects, when H is tilted 
at 8 < 8 C away from the column direction. We explain 
the critical behavior of the linear resistivity on the basis 
of the CIC transition of kink chains right above 8 C , 
as predicted in (1 + 1) dimensions for one-dimensional 
correlated disorder. This implies that the appearance 
of the TME is concomitant with the vanishing of the 
linear resistivity at 8 C . As a consequence, our results 
support the scenario for an usual BG to liquid transition 
in irradiated BSCCO in contrast to a puzzling two-stage 
BG to liquid transition recently observed in untwinned 
YBCO with columnar defects S. 
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Figure 1: V/I — I curves for tilted magnetic fields. From the 
right to the left 9 = 0, 5, 10, 15, 20, 25, 30, 35 and 40°. The 
curved lines are a fit of the Bose glass theory to the data. 



Figure 2: Angular dependence of ln(V/I) for different I. The 
solid lines are a fit of ln(V/I) = M + Ntan9 to the data; the 
dotted lines are a guide for the eye. Insert: N/I^ 3 vs / (left 
axis) and N/I\ vs / (right axis) with 7o and 7i determined 
from theory, (a) / = 1/3 and T = 68.553K (filled symbols), 
/ = 2/15 and T = 69.136^ (open symbols), (b) f = 1/3 
and T = 59.9607s: in another crystal with B$ = 0.75 T. The 
straight lines are fits whence we obtain fx = 0.33 ± 0.05 (solid 
line) for I < I2 and ji = 1.0 ± 0.1 (dashed line) for I > I2. 



Figure 3: ln(V/I) vs J- 1/3 for I < h and 9 < 9 C . Insert: 
ln{y/I) vs J -1 for I > I2 and 9 < 9 C . In both plots the lines 
are a fit to the data using the form ln(V/I) = P — Q^I -11 
with /j = 1/3 or 1 according to whether I is < J 2 or > I 2 . 
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Figure 4: Angular dependence of Q1/3 defined in the regime 
I < I2 (see Fig.^). Lower insert: angular dependence of Qi 
defined in the regime I > 1%, Upper insert: angular depen- 
dence of Qi/3/I^ 3 and Qx/I\. In each plot the straight line 
is a least- square fit. 



Figure 5: Critical behavior of the linear resistance for differ- 
ent values of / and T. The curved line is a fit to the data using 
R = R c (tan6 — tanOc) 1 ^ 2 in accordance with a CIC transition 
scenario. Insert: angular dependence of the Bose glass tran- 
sition. The solid line is a fit of T SG (if|| , 0) - Tbg{H\\ , H±) ~ 
IH^l 1 ^ 1 - with vx. = 1.0 ± 0.1 to the data, the point (0,0) 
included. 
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